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Abstract. We present a new single sum series evaluation of Moll's quartic 
integral and present two new generalizations. 

In a beautiful personal story [6] Victor Moll describes his encounter with certain 
quartic integral and derives its evaluation and goes on to study analytic and num- 
ber theoretic properties (log- concavity, p-adic valuations, location of the zeros, etc.) 
of a polynomial associated with the evaluation of the integral [1,2,4,5,6,7]. In this 
article we use the Almkvist-Zeilberger algorithm ([3,8,9,10]) to derive a new series 
evaluation of this integral. In addition, we give two new generalizations of the iden- 
tity. In [1], T. Amdeberhan and V. Moll presented a survey of old and new proofs 
of the evaluation and the formula: 

Theorem 1 [T. Amdeberhan and V. Moll, [1]]: 
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and (z) k = z(z + \)(z + 2) . . . (z + k - 1). 

The polynomial associated with the evaluation of the integral that is the subject of 
study in [1,2,3,4] is 

Next we state the main results of this article: 
Theorem 2: 
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Proof: 



We use the Almkvist-Zeilberger algorithm ([3,8,9,10]), and the reader is assumed 
to be familiar with this method. In particular, we used Zeilberger's Maple package 
EKHAD8 (procedure AZc) that computes differential operators and certificates for 
single variable hyper-exponential functions accompanying [3], available from 



http : //www .math . rutgers . edu/~zeilberg/tokhniot/EKHAD 



We cleverly construct the (certificate) function 

x(4m + 3 + Aax 2 m + 2ax 2 — x 4 ) 

RM = (x 4 + 2aX 2 + 1) 

with the motives 



-4m - 3 - 4a(2m + 3)D a (F(x, a)) - 4(a 2 - l)D 2 a F(x, a) = D x (R(x, a)F(x, a)) , 

where F(x, a) is the integrand and D a is differentiation operator with respect to the 
variable a. If we integrate both sides with respect to x on the limits of integration 
and observe that the right-hand side vanishes, we get a differential operator 

-4m - 3 - 4a(2m + 3)D a - 4(a 2 - l)D 2 a , 

that annihilates the left side of the theorem. Using the standard technique (or 
use Paul Zimmermann and Bruno Salvy's gfun from Maple library if you wish) 
of translating a differential equation satisfied by a power series into a recurrence 
relation for its coefficients a;(m), we get 



(-4Z 2 + (-8m - 8)/ - 4m - 3)oj(m) + (4Z 2 + 121 + 8)oj(m)(/ + 2) = , 

a homogeneous recurrence relation satisfied by the discrete coefficient function 
ai(m). Finally, the theorem follows by solving the recurrence relation with the ini- 
tial conditions calculated directly from the integral: ao(m) = 1(0, m) and ai(m) = 
7'(0, m), where J(a, m) is the the integral on the left. Q.E.D. 



Comparing the right-hand side of our theorem with that of (theorem 1 ), we get 



P rn (a) = 2 - ^ EC" 1 ) — , , 2 4 (Polypart) 

4-7T ^ — ' (!m! 

1=0 

Using Newton's Binomial theorem, 

(1 + ar+ l/2 = g^+l/2\ k 

k=0^ ' 

and multiplication of series, the coefficient of a™, d n (m), in the polynomial P m (a) is 
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Next, we give the first of two generalizations in which 2 in the integral of theorem 
2 is replaced by any integer n for which the integral exists. 



Theorem 3: 



f 

Jo 



rlT 1 ^2, nH l 2 "- 1 Vf ' i m 1 V 



(x 2 ™ + mx™ + l) m+1 2n ^ v l\m\ 



any integer n for which the integral exists. 



Proof: 



First, we make the change of variables z = x n and the question reduces to evaluating 

dz 
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n(z 2 + 2az+l) m + 1 z 1 - 1 / n 
Then, EKHAD gives a differential operator 

-2n - 2nm + 1 - (2m + 3)n 2 aD a - n 2 (a 2 - l)D 2 a . 
with certificate function 

, . nz(2n + 2nm — 1 + nzma + naz — az — z 2 ) 

R{z,a) = J- — 

z A + az + 1 

That is, 



(-2n - 2nm + 1 - (2m + 3)n 2 aD a - n 2 {a 2 - l)D 2 a )F{x, a) = D x (R(x, a)F{x, a)) . 

where F(x,a) is the integrand and D a is differentiation operator with respect to the 
variable a. Now integrate both sides and convert the resulting differential operator 
for the series into a recurrence relation for the coefficients and solve. Q.E.D. 

The second generalization where n is replaced by any parameter a for which the 
integral exists whose proof follows from theorem 3 by writing aa as n (^)- 

Theorem 4: 
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for any integer n > and indeterminate a for which the integral exists. 
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Problem: Find analogous polynomial Polypart as in theorem 1 associated with the 
evaluation of the generalization in theorem 3 if it exists, i.e. find dosedForm(n 7 m, a) 
such that 



P«(a) := closedForm^m^x^^y— ~ °' ' 

is a polynomial in a. 

For the special case n = 2 {Polypart), closedForm(2,m 7 a) — f (2( a +i) 1 )"'+ 1 /2 • 
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